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Abstract 

We develop the Schwinger-DeWitt technique for the covariant curvature expansion of the 
quantum effective action for brane induced gravity models in curved spacetime. This expansion 
has a part nonanalytic in DGP type scale parameter m, leading to the cutoff scale which is 
given by the geometric average of the mass of the quantum field in the bulk M and m. This 
cutoff M cuto ff = vMm is much higher than the analogous strong coupling scale of the DGP 
model treated by weak field expansion in the tree-level approximation. The lowest orders 
of this curvature expansion are calculated for the case of the scalar field in the (d + 1)- 
dimensional bulk with the brane carrying the d-dimensional kinetic term of this field. The 
ultraviolet divergences in this model are obtained for a particular case of d = 4. 

1 Introduction 

Modified theories of gravity in the form of braneworld models can in principle account for the 
phenomenon of dark energy [TJ [2] as well as for nontrivial compactifications of multi-dimensional 
string models. It becomes increasingly more obvious that one should include in such models 
the analysis of quantum effects beyond the tree- level approximation [3]. This is the only way to 
reach an ultimate conclusion on the resolution of such problems as the presence of ghosts [4] and 
low strong-coupling scale [3 El [7]- Quantum effects in brane models are also important for the 
stabilization of extra dimensions [8], fixing the cross-over scale in the Brans-Dicke modification 
of the DGP model [9] and in the recently suggested mechanism of the cosmological acceleration 
generated by the four-dimensional conformal anomaly |lUj . 

A general framework for treating quantum effective actions in brane models (or, more gen- 
erally, models with timelike and spacelike boundaries) was recently suggested in I12| , IT3" 1 I14j. 
The main peculiarity of these models is that due to quantum field fluctuations on the branes 
the field propagator is subject to generalized Neumann boundary conditions involving normal 
and tangential derivatives on the brane/boundary surfaces. This presents both technical and 
conceptual difficulties, because such boundary conditions are much harder to handle than the 
simple Dirichlet ones. The method of [13] provides a systematic reduction of the generalized 
Neumann boundary conditions to Dirichlet conditions. As a byproduct it disentangles from the 
quantum effective action the contribution of the surface modes mediating the brane-to-brane 
propagation, which play a very important role in the zero-mode localization mechanism of the 
Randall- Sundrum type [T3]. The purpose of this work is to make the next step — to extend a 
well-known Schwinger-DeWitt technique [16\ [T71 [181 [1~9] to the calculation of this contribution in 
the DGP model in a weakly curved spacetime in the form of the covariant curvature expansion. 



*e-mail: barvin@lpi.ru 
^e-mail: nesterov@lpi.ru 



1 



Briefly the method of [13] looks as follows. The action of a (free field) brane model generally 
contains the bulk and the brane parts, 

S[$] = \ J d d+1 X^G $(X)F{V x )$(X) + \ Jd&y/g y{x) x(V 2 ) <p(x) , (1) 

B b 

where the (rffl)-dimensional bulk and the (i-dimensional brane coordinates are labeled respectively 
by X = X A and x = x^, and the boundary values of bulk fields $(X) on the brane/boundary 
b = <9B are denoted by <p(x), 

0(X) = <p(x), (2) 

b 

G and g are the determinants of the bulk Gab and g^v metrics respectively. 

The kernel of the bulk Lagrangian is given by the second order differential operator F(Vx), 
whose covariant derivatives Vjc in (pQ) are integrated by parts in such a way that they form 
bilinear combinations of first order derivatives acting on two different fields. Integration by parts 
in the bulk gives nontrivial surface terms on the brane/boundary. In particular, this operation 
results in the Wronskian relation for generic test functions (fix %(X), 

J d d+1 XVG~($i~F(V x )$2-$lF(Vx)&2) = - jd d x^g ^!^ 2 -^ 2 ). (3) 

B b 

Arrows everywhere here indicate the direction of action of derivatives either on <Pi or ^2- 

The brane part of the action contains as a kernel some local operator x(V), V = V^. Its order 
in derivatives depends on the model in question. In the Randall-Sundrum model [15], for example, 
it is for certain gauges just an ultralocal multiplication operator generated by the tension term 
on the brane. In the Dvali-Gabadadze-Porrati (DGP) model [1] this is a second order operator 
induced by the brane Einstein term on the brane, x(V) ~ VV/m, where m is the DGP scale 
which is of the order of magnitude of the horizon scale, being responsible for the cosmological 
acceleration [2]. In the context of the Born-Infeld action in D-brane string theory with vector 
gauge fields, x(V) is a first-order operator [20] . 

In all these cases the variational procedure for the action ([1]) with dynamical (not fixed) fields 
on the boundary (p(x) naturally leads to generalized Neumann boundary conditions of the form 



[W(Vx) + *(V)J # b = 0, (4) 

which uniquely specify the propagator of quantum fields and, therefore, a complete Feynman 
diagrammatic technique for the system in question. The method of [13] allows one to systemati- 
cally reduce this diagrammatic technique to the one subject to the Dirichlet boundary conditions 
(P| b = 0. The main additional ingredient of this reduction procedure is the brane operator 
F hra,ne (x,x') which is constructed from the Dirichlet Green's function Gd(X, X') of the operator 
F(V) in the bulk, 

F^{x,x') = -W(Vx)G D (X,X')W(V x >) v „ v , , + x(V) x') . (5) 

X=e(x), X'=e(x') 

This expression expresses the fact that the kernel of the Dirichlet Green's function is being acted 
upon both arguments by the Wronskian operators with a subsequent restriction to the brane, 
with X = e(x) denoting the brane embedding function. 

As shown in [13], this operator determines the brane-to-brane propagation of the physical 
modes in the system with the classical action ([T| (its inverse is the brane-to-brane propagator) 
and additively contributes to its full one-loop effective action according to 

A-loop = \ Tr^ +1) InF = \ Trg +1) InF + \ Tr^ lnF brane , (6) 
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where Tr DN denotes functional traces of the bulk theory subject to Dirichlet and Neumann 
boundary conditions, respectively, while Tr" is a functional trace in the boundary (f-dimensional 
theory. The full quantum effective action of this model is obviously given by the functional 
determinant of the operator F(Vx) subject to the generalized Neumann boundary conditions 
([5]), and the above equation reduces its calculation to that of the Dirichlet boundary conditions 
plus the contribution of the brane-to-brane propagation. 

Here we apply ([6|) to a simple model of a scalar field which mimics in particular the properties 
of the brane-induced gravity models and the DGP model [lj. This is the (d+ l)-dimensional 
massive scalar field $(X) = <&(x, y) with mass M living in the curved half-space y > with the 
additional (i-dimensional kinetic term for <p(x) = ${x, 0) localized at the brane (boundary) at 

y = o, 

S[<f>] = ^ j d d+1 xVG((V x $(.X)) 2 +M 2 <P 2 (X)+P(X)<P 2 (X)^ 

+ i / ^ ^ {jy^x)? + A*V(*) + P{x)v 2 {x)) ■ (7) 

Here and in what follows we work in a Euclidean (positive-signature) spacetime. Therefore, this 
action corresponds to the following choice of F(Vj) in terms of {d+ l)-dimensional covariant 
d'Alembertian (Laplacians) 

F(V X ) = -n {d+1) + M 2 + P = -G AB V A V B + M 2 + P. (8) 

In the normal Gaussian coordinates its Wronskian operator is given by W = —d y — the normal 
derivative with respect to outward-pointing normal to the brane, and the boundary operator 
x(V) equals in terms of the (i-dimensional d'Alembertian 

x(V) = l(-D + ^+p), □ = □(<*> = g^V^u, (9) 

where the dimensional parameter m mimics the role of the DGP scale [1]. Thus, the generalized 
Neumann boundary conditions in this model involve second-order derivatives tangential to the 
brane, 

cf. ([H) with W = —d y and k given by (J9]). 

As was shown |14j . the flat space brane-to-brane operator for such a model without potential 
terms has the form of the pseudodifferential operator with the flat-space □, 

F brane (V) = — (-□ + ii 2 + 2m V / M 2 -D ). (11) 
2m 

In the massless case of the DGP model [I], M = 0, this operator is known to mediate the 
gravitational interaction on the brane, interpolating between the four-dimensional Newtonian 
law at intermediate distances and the five-dimensional law at the horizon scale ~ 1/m [6\. 

Here we generalize this construction to a curved spacetime and expand the brane-to-brane 
operator and its effective action in covariant curvature series. This is the expansion in powers of 
the bulk curvature B R, extrinsic curvature of the brane k a /s, the potential terms of the bulk P 
and brane p operators and their covariant derivatives — all taken at the location of the brane. 
The expansion starts with the approximation (jlip based on the full covariant d'Alembertian on 
the brane. We present a systematic technique of calculating curvature corrections in (J6]) and 
rewrite their nonlocal operator coefficients — functions of the covariant □ - in the form of the 
generalized (weighted) proper time representation. 
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The success of the conventional Schwinger-DeWitt method is based on the fact that the one- 
loop effective action of the operator, say — □ + M 2 , has a proper time representation 



oo 

^Trln (-D + M 2 ) = -i J ^ e~ sM2 TV e sU . (12) 

o 

In view of the well-known small time expansion for the heat kernel [16} [T7] , 

^ oo 

= -J—DV^x^e-^Vt'jr s n a n (x,x'), (13) 
( 47rs ) 1 ^ 

(a(x,x') is the geodetic world function, D(x,x') is the associated Van Vleck determinant and 
a n (x,x') are the Schwinger-DeWitt or Gilkey-Seely coefficients) the curvature expansion even- 
tually reduces to the calculation of the coincidence limits of a n (x,x') and a trivial proper time 
integration resulting in the inverse mass expansion 



-Trln [M 2 -U = — - — -— > v - j0 ' dxJga n {x,x 



(14) 



As we will show below, the calculation of the brane effective action differs from the conven- 
tional Schwinger-DeWitt case in that the proper time integral (|12[) contains in the integrand a 
certain extra weight function w(s), and instead of just Tr e sC one has to calculate the trace of the 
heat kernel acted upon by a certain local differential operator Tr (W(X7)e s ^) . This again reduces 
to the calculation of the coincidence limits — this time of the multiple covariant derivatives of 
a n (x,x'), a(x,x') and D(x,x') — the task easily doable within a conventional DeWitt recurrence 
procedure. 

The result of this calculation is peculiar. Unlike the usual Schwinger-DeWitt expansion (|14j) 
the brane effective action takes the form 

2 \^)h MN h mN ~ l 

+^£^E^-), d5) 

v I 7V=0 i<N 

where 0(9Jl k ) represent the integrals over the brane/boundary space of local invariants of di- 
mensionality k in units of mass or inverse length. With this notation, in particular, a n (x,x) = 
0(9JI 2n ). More generally, these invariants (or spacetime covariant higher-dimensional opera- 
tors) are composed of the powers of the bulk and brane curvature, extrinsic curvature of the 
brane/boundary, the potential terms of the bulk and brane operators and their covariant deriva- 
tives. 

The main difference of (|15|) from (|14p is that in addition to a usual part analytic in m with 
a typical M-dependence (second series in (|15|) ) we also have the part singular in m -> with a 
qualitatively different analytic dependence on the bulk mass 

( M d/2-N instead of M d ~ N ). This 
property was recently discovered for the effective potential in the toy model of the DGP type |14j . 
Physically this leads to an essential modification of the perturbation theory cutoff — the domain 
of validity of the local expansion 9Jt <C M cuto ff • It reduces this cutoff from M cuto gf = M to 

M cutoff = ^fM^i. (16) 

In physically interesting brane models with m M this implies essential reduction of M cuto flF 
and signifies the problem of a low strong coupling scale [6]. While in [6] this phenomenon was 
observed in the tree-level theory, here we extend it to the quantum one-loop approximation. 
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As an application of this generalized Schwinger-DeWitt expansion we calculate the one-loop 
brane effective action of the quantum scalar field with the accuracy 0(9Jt 2 ). In this approximation 
the basis of local curvature invariants includes one structure as a cosmological term, two structures 
linear in the extrinsic curvature and the potential term of the brane operator (|9|) and seven 
structures of dimensionality (Wl 2 ), 



0(2H°) = jd d x^, (17) 

b 

0(m l ) = Jd d x^k, Jd d x^^-, (18) 

b b 

0(m 2 ) = jd d x^ B R, J d d x y/g B R nn , J d d x y/g kl p , f d d x ^ k 2 , 



b b b b 

2 



b 



Here B R is a bulk scalar curvature, B R n n = b Rab n^n B is the projection of the bulk Ricci tensor on 
the normal vector n A to the brane, k a p is the extrinsic curvature of the bran^] and k = g a ^k a p is 
its trace, and P and p represent the bulk and brane potential terms of relevant operators F(Vx) 
and x(V x ) introduced above. Below we find explicit coefficients of these structures in (|15p as 
nontrivial functions of mass parameters M, m and [i and find UV divergences in this model. 



2 Perturbation theory for the bulk Green's function and brane- 
to-brane inverse propagator. 

In normal Gaussian coordinates the covariant bulk d'Alembertian decomposes as D^ +1 ^ = 
d 2 + D x (y) + where ellipses denote depending on spin terms at most linear in derivatives^ and 
^x(y) is a covariant d'Alembertian on the slice of constant coordinate y. Therefore the full bulk 
operator takes the form 

F(V) = + M 2 + P{X) 

= —d 2 — □ + M 2 — V(X | d y ,V) = F -V, □ = 0,(0), (20) 

in which all nontrivial y-dependence is isolated as a perturbation term V(X \ d y ,V) = V(y,d y ) 
- a first-order differential operator in y, proportional to the extrinsic and bulk curvatures, and 
of second order in brane derivatives V which we do not explicitly indicate here by assuming that 
they are encoded in the operator structure of V(y,d y ). In particular, it includes the difference 
n^(0) — D(y) = □ — D(y) expandable in Taylor series in y. 

The kernel of the bulk Green's function can formally be written as a y-dependent nonlocal op- 
erator acting on the cf-dimensional brane — some non-polynomial function of the brane covariant 
derivative 

G D (X,X') = G D (y,y'\V)5(x,x') . (21) 

1 In this article we use the following sign convention for extrinsic curvature 

k a /3 = e| a e^) (y x ) A n B , 

where {e^} is a holonomic basis tangent to brane, n A - unit inward normal vector to the brane. 

2 This term for a general spin structurally has the form fcVx +k 2 + (Vfc) + R where k is the extrinsic curvature 
of y — const slices and R is the bulk curvature. 
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The perturbation expansion for Gd{v, y'\ V) is usual 

oo 

G D = G° D + G° D VG° D + ... = G%J2 ( VG dT , (22) 



?1=0 



where G^ is the propagator for operator Fq obeying Dirichlet boundary conditions and the 
composition law includes the integration over the bulk coordinates, like for example in the first 
subleading term 

/•oo 

G° D VG° D (y,y')= / dy" G%(y, y")V(y" , d y ")G%(y" , y ) . (23) 
J o 

The lowest order Green's function in the half-space of the DGP model setting — the Green's 
function of Fo = —d 2 — □ + M 2 subject to Dirichlet conditions on the brane y = and at infinity 
— reads as follows 

p -\y-y'\ ^fw-u _ p ~(y+ y ') vm^n 

G° D (y, y') = ; — . (24) 

We want to stress that here we assume the exact (curved) d-dimensional d'Alembertian □ depend- 
ing on the induced metric of the brane g fiu (x). This means that in the lowest order approximation 
the underlying spacetime is not flat, but rather has a nontrivial but constant in y metric of con- 
stant y slices. Correspondingly in the zeroth order we have 

- [WG D (y,y')W] = -$,G° D (y,y')%y = ^M 2 -U . (25) 

» » y=y'=l) 

The perturbation of the bulk operator can be expanded in Taylor series in y, so that it reads 

oo ^ oo ^ 

V = k(x, y)d y + U(y) -a- P(x, y) = ^ -u k y k d y -J2 y k , (26) 

fc=0 ' fc=0 

where itfe(V) and ffc(V) form a set of y-independent local d- dimensional covariant operators of 
maximum second order in V x - The coefficients of these operators are given by the powers of the 
bulk and brane curvature, the extrinsic curvature of the brane, the potential term P and the 
covariant derivatives of all these quantities — all of them taken at the brane. 

Below we present them up to the first order in the bulk and brane curvature and to the second 
order in the extrinsic curvature of the brane. Working in gauss normal coordinate system for the 
case of a single quantum scalar field they read as 

u = k(x,y)\ y=Q = k, 

«i = d y k(x, y) | y=Q = - B R nn - kl„ (27) 

and 



y=0 



v = P{x,y) 

v 1 =d y (- U x {y) + P(x, y)) = 2k a P V a Vp + 2(V Q fc a/J ) - (V^fc) + 0[M% 

y=0 

v 2 = d 2 y {- D x (y) + P(x,y)) | = (-2 B R a f n -6k^k/) V a V p + Q[Tl 3 ]. (28) 



Here we everywhere omit the argument x of all quantities located on the brane, V a is the d- 
dimensional covariant derivative on the brane, k a n denotes the extrinsic curvature of the brane, 
k = g a ^k a /3 is its trace. Subscript n denotes the projection of the relevant bulk index to the 
normal vector n A , in particular B R nn = B R c ACB n A n B . 
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Another source of differential operators with coefficients of growing power in the curvatures 
and their derivatives is the commutation of □ with all x-dependent quantities involved, like 

[ □, no ] = P, k(x)\ = 2 {V a k)X7 a + 0[9Jl 3 ] 

[D,vi] = 4(V^)V^V a V /3 + 0[9Jt 3 ]. (29) 

In all these equations 0[9Jt'] denotes the accuracy in powers of the dimensionful quantities 
— curvatures, operator potential terms and their derivatives — with which the relevant quantity 
is calculated within the local Schwinger-DeWitt technique. In fact I — 1 associated with 0[VJl l ] 
above indicates the dimensionality (in mass units) of the coefficient of the higher-derivative term 
of the relevant operator — its symbol. In [17J it was called a background dimensionality, in 
contrast to a total dimensionality of the quantity. In what follows we will denote the background 
dimensionality of the operator O = 0[9Jt l ] by 

DimOpSm*] = 1. (30) 

Thus for an operator of the form (V a R m P l k n )V b , where the coefficient of the 6-th order derivative 
is given by a monomial of curvatures and potential terms and their covariant derivatives of the 
total order a, the background dimensionality is given by 

Dim (V a R m P l k n )V b = Dim o(M a+2m+2l+n ) = a + 2m + 2l + n, (31) 

whereas its total dimensionality is, of course, [(^/ a R m P l k n )^J b ] = a + b + 2m + 21 + n. 

The background dimensionality in fact counts the order of the perturbation theory in powers 
of the local quantities R = 0(W1 2 ), P = 0(T1 2 ), k = Offi 1 ) and their derivatives. The orders of 
this perturbation theory have the form 9Jt N /M , where M is the bulk mass playing the role of the 
cutoff, fJJV^M, beyond which the local expansion does not apply. The background dimensionalities 
of relevant local operators, like (f28|) - ([29|) . always contribute to N. On the contrary, free (acting 
to the right) derivatives of these operators may effect the overall power of the cutoff M in the 
denominator (by reducing it) so that they are not indicative of the accuracy of the Schwinger- 
DeWitt expansion. In other words, any quantity of the background dimensionality I contributes 
to the Z-th order of the local expansion 9J1 1 and higher, whereas the total dimensionality of this 
quantity does not determine the order of this expansion. This explains a distinguished role of the 
background dimensionality vs the total one. 

The next calculational step consists in the substitution of (|25p and (|26p into (|22p , and it leads 
to exactly calculable integrals over y. The integration over y results in a nonlocal series in inverse 
powers of \/M 2 — □ — this is obvious from the y-expansion of (|26p . because every extra power of 
y brings one extra inverse power of v 7 M 2 — □ . Each k-th order of this series arises in the form of 
the following nonlocal chain of square root "propagators" , 

— ; 1 — V\ — ; 1 — V 2 ...V v -i — ; 1 , , h + h + - + h = k, 
(VM 2 -^ (VM 2 -D)fe P (^M 2 -Oyv P 

with some differential operators Vi as its vertices. With the aid of the commutation relations like 
(|29|) all these propagators can be systematically commuted to the right of the expression by the 
price of extra commutator terms of the same structure, and the perturbation expansion finally 
takes the form 

oo 1 

- [WG D (y,y')W] y=yl=Q = - g^V)^—^, (32) 

where £/fc(V) is a set of certain local covariant differential operators acting on the brane. The 
dimensionality of each L^(V) is the inverse length to the power k, which is composed of the 
dimensionalities of bulk and extrinsic curvatures and covariant derivatives all taken on the brane 
at y = 0. 
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With x(V) given by ([9]) the brane-to-brane operator ([5]) reads 



(V) = F -E^(V) (VM2 ' D)fc _ i ; (33) 

F = f - □ + // 2 + 2m x/M 2 -^ ) , (34) 

Here we absorbed the potential term p of the operator x(V) into the first term, k = 1, of the 
perturbation part of _F brane (V) by redefining the U\ term, C7i — ► U\ — p/2m, because p should of 
course be treated on equal footing with other perturbation:^ 3 ]. The rest of the perturbation part 
is induced from the bulk and does not depend on the boundary conditions on the brane encoded 
in the operator x(V). Note that the zeroth-order term here is a nontrivial nonlocal operator 
because □ is a curved space (i-dimensional D'Alembertian acting on the brane. 

In the 9JT 2 approximation (involving the terms linear in the bulk curvature and potential P 
and the terms quadratic in the extrinsic curvature k and the brane potential p) the operator 
coefficients t/fc(V) extend to k = 6. Higher order coefficients go beyond the 9Jt 2 -approximation. 
The calculations show that for a single scalar field they read as 

y ' 2 u 2m 2 2m 
U 2 {V) = —g ui - -v - - Uq n = - B R nn + g fcj, - - P - ^ k 2 , 

= ~\ k al3 V a Vp - X - k^.aVp + 0(W1 3 ), 
U,(V) = ~\v2 + 0(9Jt 3 ) = (1 B R\? n + \ k^k/j V a V p + 0(9Jt 3 ), 

Us{V) = - 1 [□ , Vl ) + o(m 3 ) = -\{V»k a P) V M V Q V^ + O(0JT 3 ), 

o Z 

Ue(V) = Jr i>i «i + 0(W1 3 ) = ^ k^k^VaVpV^Vv + 0(W1 3 ). (35) 

As mentioned above, each Uk has a total dimensionality k in units of mass. Except the case of 
k = 2, for which U2 = O(50T 2 ), their background dimensionality is Uk = 0(9Jt^ fc+2 - > / 3 ^), where the 
square brackets denote the integer part of a fractional number. 



3 Perturbation theory for the brane effective action 

Perturbation theory for the effective action immediately follows from the perturbation series 
for the operator F brane . The brane effective action can be rewritten as 



- TrlnF branc = -TrlnF + -Trln ( 1 - V U k (V) , = - ) , 



(36) 



and reexpanded in powers of the perturbation series term under the logarithm sign. After com- 
muting all the square root "propagators" 1/(V M 2 — □ ) k and the propagators 1/Fq to the right 
this expansion takes the form 

-TrlnF brane = -TrlnF -- V Tr^(V) ^ 3_, (37) 

2 2 2 ^ \VW^UfF^ 



3 We do not introduce a new notation for Ui, and this should not lead to a confusion because Eq. (|32|l will not 
be used in what follows. 
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with a new set of local covariant differential operators Wm(V) acting on the brane. For dimen- 
sional reasons the total dimensionality of Wj~i(V) is k+l in units of mass, because Fq like \JM 2 — □ 
has a unit dimensionality. 

These operators are composed of the products of the operators C/fc(V) introduced above and 
their multiple commutators with the the "propag ators" 1/(VM 2 -D) fc and 1/F . These commu- 
tators are based on the multiple use of the formula 



/(□)£ - Bf(D) = E fci ad n B ■ a 5/( D )' ad D 5 = P' B l 

k=l 



(38) 



which leads to the following structure of W^iCV) 



max{0, 1-2} 
p=0 



(39) 



Here, modulo the powers of p/2m originating from C/i(V), the coefficients Wfcz iP (V) are m- 
independent, and the negative powers of m follow from the differentiation of the propagator \/Fq 
with respect to □ participating in the commutators [Uk, l/-Po], 



l 



1 1 
2iV 



1 1 



-□ 



As we see, each such differentiation results in the extra power of the propagator 1/Fq and an 
extra term proportional to 1/m. For a Wki(V) with a given / the highest order of 1/m is limited 
by max{0, 1—2}. This is because p differentiations increase the power of 1/Fq from some initial /' 
to at least / = I' +p, and the initial /' > 2, because the commutation of 1/Fq with other quantities 
begins only when the number of these propagators exceeds two (one should remember that one 
propagator 1/Fq always stands to the right of everything else, see Eq . ([36]) ) . This explains the 
upper limit of summation over p in (|39p and, as we will later see, underlies the regularity of the 
Neumann limit m — > oo. 

The background dimensionality of W^^ P (V) is a monotonically growing function of all its 
indices and can be shown to satisfy the bound 



DimW W)P (V) > 



k + I + 2p + 2 



(40) 



where square brackets denote an integer part of the fractional number. This bound follows from 
the observation that for any Wki, p , composed of the chain of and n commutators with □, 
the total dimensionality = k+l+p equals [Wm jP ] = Yli ^ + ^n. On the other hand, the 

background dimensionality in addition to the sum of DimC/^ = [(ki + 2)/3] contains at least one 
extra unit of mass per each commutation with □. Therefore 



DimWiH,p(V) >Y, 



'h+2' 


+ n > 






+ n = 


3 




3 



{k+l+p)+n+2 



(41) 



where we used the above counting of the total dimensionality ^ ki + 2n = k+l+p. The bound 
(|4U|) then follows from the fact that the overall negative power of m does not exceed the number 
of commutations n, < p < n. 

Relevant Wki which can contribute up to 0(pJt 2 ) order inclusive are 



Wi 



01 



W 



ii 



Uo 



-k 

2 



P 

2m' 
1 



P 
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2 

^ = -!^ = - t-k + ^ ' 

2 1 2\2 2m, 
W 2 i = U 3 = -\ k afi V a Vfs - i {V a h^)Vp + 0(9Jt 3 ), 

w 3 i = u 4 = (~ fl i? Q / n + 1 v a v p + o(9Jt 3 ), 
w 22 = u 3 U! + \u 2 u 2 = \(h + JL\ k a ? v a v p + o(m 3 ), 

2 2 \2 2m J 

w 4 i = u 5 = ~{Vk a ^ v^VaVp + o(m 3 ), 

W51 = U 6 = l k^k^VaVpV^v + 0(W1 3 ), 
o 

W 42 = UsXJx + U A U 2 + ~U 3 U 3 + U 3 [D, Ux] = i k^krVaVpVv&v + 0(^ 3 ). (42) 
4 Generalized proper time method. 

The further calculation is based on the possibility to express the nonlocal structures in (|37p in 
terms of the heat kernel of the box operator □, which admits a well-known curvature expansion. 
This is the set of proper time representations which differ from a usual Schwinger integral by 
nontrivial weight functions wm(s) 

00 

InFo = -ln(2m) - / ^-w 00 (s) e"^ 2 " ) , (43) 





00 

I I f ds -s(M 2 -0) 



-7 = / -w kl (s)e-^ M ~^>, (1>1). (44) 



(y/M 2 -D) k F l J s 


These weight functions can be found as follows. First, decompose the operator Fq defined by 





i a 

(|34p into the product of two factors linear in V M 2 — □ 

^-(VM 2 -n - m + )(y / M 2 -D - m_) , (45) 

where m± denote the roots of the quadratic equation x 2 + 2ms — M 2 + fi 2 = 0, 

m + = —m + \J m 2 + M 2 — m_ = — m — \J m 2 + M 2 — fi 2 . (46) 

Therefore 



F 

2m 



In F = - In 2m + In (V 'M 2 -U - m ) + In {\fM 2 -U 



m 



(VM 2 -D) k F Q l ^ WM 2 -ny 



I 



E( Dl, (to, , m 1 1 



^ \(\/M 2 -n -m V- 
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Here the second equation is the result of the decomposition of its left hand side into partial 

d 

(-m/2) 1 1 



fractions with the coefficients B?, and 



x 



' r(2l-b) r(k-a+b) f 1 ( ' m^m_\ b 1 fm-m^ 1 ' 

m k ~ u \ m 



and 



Dj;,(m m ) = (™V -1 t£. g A*+'-M ^ ) . (49) 



A / r(k) m + k+l ~ b ^ r{l-b+l-p) p\ \m-m_ 
Now, in addition to a simple proper time representation of the square root "propagators" 

CO 

1 f ■ a/2-1 -s(M^-U) 



(VM 2 -D) a r(o/2) 



da^/^e- 8 ^ - u> , (50) 



we need a similar representation for the new operator of the form \J M 2 — □ — m and its logarithm. 
It can be derived with the aid of the integral representation of the cylindrical function of a 
half-integer order 



OO 



(VM 2 -D-m) a Pi 



a J 



o 
o 



20Fr(l + a) 

u 

Differentiating it with respect to a at a = one gets 



s " 3 / 2 e - s (^ 2 -°) J dx x a e mx e~£. (51) 
o 



CO CO 

2 



m(VM^O-m) = -^j= J ds s~ 3 / 2 e"^ 2 " ) J dx e mx e"t 



CO 

1 / ' d £ u; (_ m ^) e -(M 2 -n) ! (52) 



2 

where the function w{—m^fs) is given in terms of the complementary error function erfc(z), 

CO 

2 



w\—a 



L f dx e 2ax e~ x2 = e a2 erfc (-a), (53) 
o 

CO 

2 f 2 

erfc(z) = -= dte' 1 . (54) 

z 

A multiple differentiation of (|52[) with respect to m then gives 

e -( M2 - n ). ( 55 ) 



CO 

2 



1 _ 1 [ dss a/2-i d a w(-a) 

(VM 2 -a-m) a 2r(a) J da a 



a=mWs 
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Using (|50p . (j52[) and (|55p in (|47|) we finally come to the following expressions for the weights in 
Eqs.fPH 



woo = 2 ( w {- m +Vs) + w(-m_yfs) ) . 



B^(m + ,m_) /2 



2I t a / 2 ) 

In terms of these weights the action (j37|) takes the form 



+ {m + <-» m_} . 



oo 

1 Tr lnF b r ane _ _1 £ jf ^ ^ Tf ^ > 



M=0 o 



(56) 



(57) 



(58) 



where we disregarded the contribution of the local measure -Tr ln(2m) ~ <S<#(0) and, in addition 
to W M of d37]), introduced 



Woo = 1, W k0 = 0, k>l. 



(59) 



5 Generalized Schwinger-DeWitt expansion. 

The calculation of (|58p is based on the heat kernel expansion (|13p for the covariant d'Alember- 
tian □ acting in a curved <i-dimensional space without boundaries. The Schwinger-DeWitt coeffi- 
cients a n (x, x') in this expansion represent brane curvature invariants of the growing background 
dimensionality 2n, a n (x,x') = 0(9Jt 2n ). With this expansion the functional traces in (|58p take 
the form 



Tr(w kl {V)e sU ) = ( d/2 ^2 s n d d x ^ W kl (V(s)) a n (x, x>) 
nS n=0 b 



where the new generalized covariant derivative 

V a (s) = V Q - ld>(x,x') + ld x a lnD(x,x') 
Is I 

originates from the commutation of V a with the exponential factor exp(— a(x, x')/2s) and the 
Van Vleck-Morette determinant D{x,x') in the kernel of exp(sD). In addition to 

V m ...V M a n (x,x) 

x'=x 

this lengthening of the covariant derivatives also generates the coincidence limits 



V Mi-" V Mp Cr ( X > 3;/ ' 



• V m ...V Mp L>(x,x') 



easily calculable by the DeWitt recurrence procedure [HI [T7]. Moreover, this leads to extra 
negative powers of the proper time, so that 



d d x^W kl (V x ) a n (x,x') 



max{0,/— 2} [ g J 

E E ' 

p=0 c =0 



(A 



n)kl,p 



s c mP 



(60) 
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where {A n ) c klp represents the set of integrals of bulk and brane curvature invariants, powers 
of potential terms P and p/m and their derivatives of the growing dimensionality k+l+p + 
2n— 2c (which now coincides with their background dimensionality, because they are no longer 
the differential operators), 

(A n ) c kl ^ = 0(Tl k+lWn - 2c ). (61) 

The highest power of 1/s in (|60p is determined by one half of the order of W^ iP (V) in derivatives 
and equals the difference between the total dimensionality of this operator and its background 
dimensionality. As W^i is a sum over the powers of the inverse DGP mass scale (|39p . the quantities 
(An)%i p above represent the relevant coefficients of l/m p . 
Thus 

l TrlnF bra ne= 1 _}■ y [ds d/2 () -.M-C^OSy, ( g 2) 

{fc,i,n,p, c) g 

where the domain of summation over all indices is given by 

oo max{0, 1-2} [ g J 

E = E E E (^) 

{fc, Z,n,p, c} k,l,n=0 p=0 c=0 

The quantities (.An)£z p play the role of integrated generalized Schwinger-DeWitt coefficients. 
For a scalar field in the approximation 0(9Jt 2 ) only the following p = coefficients contribute to 
the brane effective action 



C4o)oo,o = Jd d x^, 

b 

G4l)8o,o =J<&y/g l h R(g) = fd d x^ {^ B R-\ 

b b 

(A)8,„ =/AvS (-^-i)' 

b 
b 

(Aft, ^/AV5 + + 

b 

(^0)21,0 = y<^vs ^ 

b 
b 
b 

G4>)§i,o = jd^Vg (J^ k lfs + ^ k 



Rnn ~\~ — k — k, 



a/3 j > 



(64) 
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6 Large mass expansion and its cutoff scales 

Integration over s in f)62[) gives 



l Trlrl pbrane _ 1 M d ^ c (Ajfc^g 

y ' {k,l,n,p,c} 

where Ci, (with j = n — c) are the following functions of the mass parameters of the model 

oo 

Cl = M 2 i- d+k+l I ^ S i- d ' 2 w M (s) e~ sM2 . (66) 
o 

The behavior of these functions for M — > oo is important for the determination of the efficiency 
of the expansion ([55]) and of the range of its validity — the cutoff M cuto g below which, Tl^M cuto s , 
this expansion makes sense. This behavior easily follows from a simple observation that the 
functions Ci, can be directly obtained from the nonlocal form-factors (|44p by integration over 
their argument □ = —A with the weight \ d l 2 ~i~ l . In the domain of convergence of this integral 
in the complex plane of d we have 



oo 

Q yd/2-j 







A (VM 2 -D) fc F l 



^-^CL. (67) 



Hf2j-d+k+l kl 
□=-A 



With the replacement of the integration variable x = \/l + X/M 2 — 1 the integral representation 
for Cl, takes the form 

o 

<p[x) = { —L> (68) 

where e± = 1 — m±/M. For M — > oo the parameter e + — > (e_ — ► 1), and the integral here has 
a nonanalytic in e + part because of the singularity of its integrand at x = [14 



darz"" 1 ^ + e)"' p(z) = e ^^f V \ (jS) + 0(1) . (69) 

1 [I) 
o 

Since e+ — * m/M in this limit, we have 

Thus, in view of the background dimensionality of {A n ) c kl = 0(Wl k+l+p+2n ~ 2c ) the local 
expansion (|65p for the effective action takes the form 

1 TrlnFbra nc = fMrnY' 2 0(9?I fc+ ' +p+2 "- 2c ) 



2 I 47T / ^ J^n-c+fc+i 

{£:,£, n,p, c} 



M d m l-p (f m k+l+p+2n-2c} 

+ (AttYI 2 J^2n-2c+k+2l ' 

V ; {k,l,n,p,c} 
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By introducing the summation index N = n — c+k + l — an overall power of 1/M — in the first 
sum and correspondingly the summation index L = 2n — 2c + k + 21 in the second sum, one can 
rewrite this series in the form (|15p presented in Introduction 

d/2 oo , n(mt^N-i\ 



1 Trln /^ane = ( MmY 1 ^ ggg 

2 " I 4tt y 2^ M N 2^ m N-, 

N=0 i<N 



M d 



oo ^ 

ExfrE^°(^" 1 )' ( 72 ) 



where the coefficient of any power of M turns out to be a finite sum of terms of a limited order 
in 5DT — the background dimensionality of relevant field invariants. 

The finiteness of such sums over i = k + I — p in the first series of (|72[) follows from the 
following simple argumentation. The range of summation over c < [ (2k + 21 +p)/6] in (|7ip is 
less than [ (2/c + 3/)/6] because p < max{0, 1 — 2}. Therefore, 



k + l = N + c-n<N + c<N + 



k l__l 

3 + 2 ~ 3 



<N + ^(k + l), 



so that k + / < N , and the ranges of summation over c, p and i = k + 1 — p indeed turn out to be 
limited at least by N. 

Similarly, the finiteness of sums over i = I — p in the second series of ([72]) is based on the 
following chain of inequalities 



k + 2l = L + 2c-2n<L + 2c<L 



2k , 2 

h I 

3 3 



(73) 



so that [k/3 + I] < L, and the ranges of summation over c, p and i = I — p are again restricted 
from above for any given L. 

This property is very important for the efficiency of the perturbation theory with the cutoff 
scale M, because otherwise any given order in 1/M would require an infinite series in rn 
- the price one could have paid for the presence of the second scale m. Fortunately, for any 
N only a finite order 0(9Jt 2N ) of perturbation theory is required. This follows from a special 
asymptotic behavior of the coefficient (|70p which brings extra powers of 1/M to (|65p . The form of 
the asymptotics (|7U|) is responsible for the two series in the expansion (|72|) . having qualitatively 
different analytic behavior in M and m, — the property recently discovered for the effective 
potential of the toy DGP model [13] • Whereas the second part is analytic in a small DGP scale 
m — > 0, the first "nonanalytic" part is formally singular in this limit, and this leads to the 
redefinition of the cutoff M cuto fj of the theory, below which 971 <C M cu t ff the local expansion 
remains valid. 

Indeed, despite the efficiency of the obtained expansion, in the first series of (|72p it contains 
negative powers of the DGP scale m and blows up for small m — > 0. This is a typical situation of 
the presence of a strong-coupling scale [6j. In fact, for m < M the actual cutoff is lower than M 
and is given by the expression (|16p . M cu t ff = V Mm, presented in Introduction (the condition of 
smallness of the strongest i = term in the first series of (|72p . 9Jt 2 /Mm ^1). 

The actual calculation of the functions (|66p can be done by using the proper time weights 
(|56p and (|57p . Since these weights imply explicit symmetrization with respect to m±, they take 
the form 

C L = ci l (M,m + ,m-) + (m + ^ m_). (74) 
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In the 0(9Jt 2 )-approximation the relevant C 3 kl turn out to be 



^01 



III 



(m 



-)$ {0) (2j + l, 2), 



m\ 2 



-_d> (o) (2j + l, 2)+$ (0) (2i + 2, 3)), 



CL = 2Q* w (23 + 2, 1), 



V A/ (2)V 



* (1) (2j + 2, l)-$ (2) (2i + 2, 1) 
+ $ (i) (2j + 3, 2)), 

Cg'i = 2(?)*(^' + 2 ' ^ 



= 2 *(")( 2 ' + 2 ' X ) " 3 *C«)( 2 J + 2 ' X ) 



V A 



2M 



A 



+ $ (3) (2j + 3, 2)) 



Here 



A= 



m.—m 



the basic function <fr (0) (a, 6) is given by the regularized Gauss hypergeometric function 
$(o)( a ' & ) = $(o)( a > b \ a ) /9M 

<T=m + /2M 

a+6+l-d 



= r( a - d) r(b) 



-; <r + 



(75) 



(76) 



(77) 



a=m + /2M 



and 



1 n_1 1 



fc=0 



a 



cr=0 



(78) 



a=m,/2M 



is the function $ (0) (a, 6|cr)/cr n with the singular at a = part subtracted, also taken at a = 
m + /2M. 

The transformation property of the hypergeometric function from the argument z to 1 — z 
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e+ \ — 2 — ^ a + o — 1 — a 



allows one to rewrite (|77p in the form which reveals the structure of the expansion ([72 

d+l-a-b 

+ 



( o)V > / v 2 / V 2 

d + 1 — a — 6 a + b + 1 — d d + 3 — a — 6 e + 
2 ' 2 ' 2 ' ~2 



r( a - d) r(b) r(d=^±r 



^fd + l — a — b a + b+l — d d + 3 — a — b e+\ 

2 ' 2 ; 2 ; fj' (79) 

where e+ = 1 — m + /M as in (|70p. In view of e+ ~ m/M the first term here generates the first 
nonanalytic in m series of (j72|) and the second term is responsible for the analytic part because 
the hypergeometric function is expandable in Taylor series in e+/2 — > 0_ 

6.1 The brane effective action in the 0(9Jt 2 ) approximation and UV diver- 
gences 

Substituting the curvature invariants of (|64p into (|65|) we get the lowest orders of the brane 
effective action in terms of the curvature invariants (|17j) - (|19j) listed in Introduction, 



1 M d I' 

l Trlnjpbr ane = _____ J ^ ^ q 

b 







b 



' 2 (4vr) d / 2 

b 

1 M d ~ 2 
2( _ J - 



+ ^~ g 0)0 + 7 ^11 ~~ g ^31^ B Rnn 

+ (~6 C ° lo+ 4 8 16 C5 ~ l2+ 16 C ^ ^ 

+ G c ° lo ~ \ c ° u+ 1 c ° o2 ~ I c ™ + h C ^ 2+ h k2 

+ G^4^^ + ^^) 2 } + °H- (80) 

Here the coefficient functions C J kl are given by equations (f~|) - (f~|) and represent a set of very 
complicated functions of M, m and /x. One can check that for [i = the Cq term given by 
Cqq = $ (0) (0, l) + (m + <-> m_) coincides with the effective potential calculated for a toy DGP 
model in [14J. 

More instructive are the ultraviolet divergences of the brane action which we present here for 



4 The a = m_/2M part of the action originating from the second term of (|74|l contributes only to the analytic 
part, because m_/2M — » — 1/2 and the relevant large M expansion originates directly from the representation (|77|l 
which does not give rise to nonanalytic terms. 



17 



the four-dimensional case in the dimensional regularization d — ► 4. They read 

1 XrlnjF brane|div = _J__^ J dXy fg (_ 4 m 2 (M 2 + 2m, 2 - 2^) ~ ^ 

b 

+ 32^ 2 (4-rf) / Q m ( M ' + 12m ' " 3 ^ 2 ) k + " 4m2 )^) 

b 
b 

+ (^ 2 -^ 2 )^+(-2m 2 + i / , 2 )fc 2 
-4m 2 P -^mkp-p 2 ^j + 0(W1 3 ). (81) 

This result confirms the general properties of ultraviolet divergences in any dimension d. 
These divergences are contained in both series of the expansion (|72|) . For an even d they are 
analytic and polynomial in both M and m, for an odd d they have a structure V Mm times 
a finite polynomial in M and m. Finally, their background dimensionality is always bounded 
by 0(9Jt d ). These general properties follow from the property of the integral (|68p which is UV 
divergent at the upper limit only for d > 2j + 21 + k. Therefore, the background dimensionality of 
the relevant terms in (fTTj) (with j = n — c) satisfies the bound k + l+ p + 2n — 2c < d + p — I < d, 
because p < max{/ — 2, 0}. The relevant overall powers of m and M in the nonanalytic part of 
(|7ip are also positive, because for the same reasons p+n — c < d/2 and n — c+k + l < d/2. Finally, 
in the analytic part of (|7ip the overall power of M in the divergent terms is again nonnegative, 
because 2n — 2c + k + 21 < d. 



7 The Neumann and Dirichlet limits 

As we see, the curvature expansion in brane induced gravity models is essentially more com- 
plicated than for pure Dirichlet and Neumann (Robin) boundary conditions. Even the conformity 
of our results with these two limiting cases (corresponding respectively to m — > and m — > oo 
in ()10|) ) requires nontrivial calculations. Here we present these calculations and check the consis- 
tency of the Neumann limit to the 0(5DT 2 ) order in the curvature, and verify the Dirichlet limit 
to all orders in 9Jt. 

First we present the known results for a local inverse mass expansion for pure Dirichlet, 
<I>|b = 0, and Robin, (d n — S')$|b = 0, boundary conditions (see [TO] and references therein). In 
these two cases labeled respectively by D and N this expansion for the effective action in the 
(d + l)-dimensional bulk reads 

oo 

- Tr^ InF = -- [ i_ e~ sM ' 2 V s n / 2 A°/ N 

2 °I N 2 J s (4vrs)( d + 1 )/ 2 ^ 

o n ~° 

1 M d °° r( ~d~l+n \ 

- 1 M V jV 2 ) A D/N / 82 n 

2 (4vr)( rf + 1 )/ 2 ^ M n_1 n y ' 

y ' n=0 

where A® represent the bulk and boundary integrals of the relevant Schwinger-DeWitt coeffi- 
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cients. The first four of them for the Dirichlet case read 



Ag = J d d+1 xVG, 

B 

b 

A 2 = j d d+1 xVG (-P+l»Rj + J d d x^g (~k\ , 

B b 

^ =^jd d x^g{+\p-^R + ^ B R nn -^e + ^kl p ). (83) 



5_ 

192 '" ' 96 

For the Neumann (Robin) case together with the bulk curvature and the extrinsic curvature 
of the boundary they involve the coefficient function S from the Robin boundary condition, 



A" = j d d+1 X VG , 

B 



A$ = J d d+1 xVG (-P+l B Rj + J d d x ^g (-~k-2S 

B b 



+ W2" + V^ + 2 kS + S ) ■ < 84) 

7.1 The Neumann (Robin) limit. 

The Robin limit of the boundary condition (|10|) corresponds to m -> +oo with M/m = 0, 
///m = and the finite limiting value of p/2m — > 5. This implies the following limits for the 
auxiliary mass parameters 

m 

m + — > , m_ ^ — oo , A — > +00 , — — > 1 . (85) 

Therefore the contribution of <r = m^./2M — > — oo terms in the coefficients C J kl (M,m ± ) given by 
(|74|) -(fT8|) vanishes in virtue of the asymptotic behavior of the hypergeometric function J? 1 (a, b; c; z) 
at z — ► — oo. Moreover, only the functions .(a, 6) with a special combination of the hypergeo- 
metric function's indices (|77p arise, and due to the relation 

they equal 



2 2 2' 2 y r(§+i)r(|+i) 



' j n! r(|+§+i)r(|+2+i) n\r{\) l2+2j l2+2j ' 

where we used the gamma function identity r(2z) r(l/2) = 2 2z ~ l r(z) r[z+l/2). As a result 
the coefficient functions CL reduce to 



k+l-d\ 



C 3 oo = ^nJ-d/2), Cl = ^ 2 ' ,fc + *>0, 
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and the expansion (J80J) for the case of the Neumann boundary conditions takes the form 

i™*— = 4 " £iL£*U, ( 87) 

v ' n=0 

where the first three coefficients equal 

Ao = ~ jd d x^ , 

b 



b 

2 4 V2m7 2 \2m) 

The consistency of this result with the Robin case is based on the duality relation ([6|) which 
implies the following identities for the coefficients of the Schwinger-DeWitt expansion (|82p 

A = -/LK +1 -^ +1 ). (89) 



These identities can be directly checked for ([8"5|) -(|84j ) with S = p/2m for n = 0, 1, 2. 



7.2 The Dirichlet case 

The Dirichlet limit formally corresponds to m — * when m/M — > 0, m//i — * and 

m^iVM 2 -^. (90) 

This limit cannot however be directly taken in the expansion (|72|) because of its obvious non- 
analyticity at m = 0. The reason is that for small m the behavior (|70p of C?, used in the original 
curvature expansion series (|65|) no longer applies. Indeed, when mM < /x 2 the parameter e + in 
the integral (|68p has another asymptotics /i 2 /2M 2 independent of m, and C£, ~ m for m — > 0. 
Therefore, in this range of m the DGP scale arises in (|65p only in positive powers I — p > 0, 
because p < I — 2 — the curvature expansion takes the form qualitatively different from (|72|) . 
The only nonvanishing term in the sum over k and I is the one with k = I = (cf. equation (|59p ). 
The relevant coefficients (|75p are 

C^(M,m ± ) = \J ^( W (-mV5) +w(-mVi)) ^ = , (91) 



because in the Dirichlet limit m 2 = m 2 = M 2 —/j, 2 and the error functions in the definition of 
the proper time weight w(— ijia/s) satisfy erfc (— m ± y / i) = 1 db erf (y M 2 —\j? y/s). These are the 
prefactors of the integrated Schwinger-DeWitt coefficients A2j = f d d x aj(x,x) in the expansion 
of the d-dimensional effective action for the massive operator ([9]) 

i TV in *(V) = i TV In (-□ V+p) = -\ jjfa g ifi^ . (92) 

This can be seen already at the level of ^-expansion (|36p . In the limit m — * all perturbation 
terms in (|36p with 1/Fq = 2m/(^ 2 — □) vanish except the contribution of U\/Fq = —p/{n 2 — □), 
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which complements the first term of this expansion to theone-loop action of the full brane operator 
with a potential (|92p (we of course disregard the volume divergent Tr In 2m part canceled by the 
local measure). 

Thus under the m — > limit the Neumann-to-Dirichlet reduction ([6]) leads to 

lim - Tr { ^ +1) In F = - Tr ( „ +1) In F + -Tr In x . (93) 
rrwO 2 N 2 D 2 v ' 

The second term here originates from the definition of the generalized Neumann vs Dirichlet 
functional determinants of the operator F 

exp(-iTr^ +1) lnF) = J D$ exp ( - S B [$] ~ \ Jd d x y x(V) , (94) 

b 

exp(- l - Trg +1) lnF) = J D$ exp ( - S B [$}) S{ <p(x) ) , (95) 

where <Sb[^] is the bulk part of the action ([1]) and the expression for the delta function of 
(p(x) = $(X) \ B (regularized by m -> 0) reads as 

5( <p{x) ) = lim Q (Det xj ^ exp ( - i J d d x y/g ^x(V)^) . (96) 

b 

The preexponential factor here, (Det xj = exp (Tr lnx/2), gives rise to the last term in ([93]) 
and, thus, confirms the consistency of the Dirichlet limit. 



8 Conclusions 

Thus we have constructed the covariant curvature expansion in massive brane induced gravity 
models, found its peculiar structure (|15j) nonanalytic in the DGP scale and derived a nontrivial 
cutoff (|16p of this general expansion. Finally, we calculated several lowest orders of this expansion 
for a quantum scalar field in a curved bulk spacetime with a kinetic term on the brane to a 
quadratic order in background dimensionality and found its ultraviolet divergences for the case 
of a 4-dimensional brane. 

These results might find important applications. Although a comparison of our massive model 
with the massless DGP model of pQ is not straightforward, we can observe a common feature 
in their cutoff properties. In both theories their cutoff (|16p is different from the bulk one M 
and is modified by the DGP scale m. For the tree-level DGP model with the Planck mass 
M = Mp, playing the role of the bulk cutoff, this cutoff equals M cuto fj = (tt^M) 1 / 3 [6j. With m 
identified with the cosmological horizon scale, this is about (1000km) -1 which is much below the 
submillimeter scale capable of featuring the infrared modifications of the Einstein theory [21J. As 
we see, the situation with the local expansion for the quantum action is much better — the cutoff 
()16|) is a geometric average of M and m, which is much higher, 

(mM) 1/2 » (m 2 M) 1/3 , (97) 

and comprises (0.1mm) -1 . This supports the conjecture that the replacement of the weak field 
perturbation theory by a derivative expansion [22], as is the case of the local Schwinger-DeWitt 
series, might improve the range of validity of the calculational scheme. 

Obviously, the Schwinger-DeWitt technique in brane induced gravity models turns out to be 
much more complicated than in models without spacetime boundaries or in case of boundaries 
with local Dirichlet and Neumann boundary conditions. It does not reduce to a simple book- 
keeping of local surface terms like the one reviewed in [19]. Nevertheless it looks complete and 
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self-contained, because it provides in a systematic way a manifestly covariant calculational pro- 
cedure for a wide class of boundary conditions including tangential derivatives (in fact of any 
order). On the other hand, the calculational strategy of the above type requires a further ex- 
tension, because there is still a large set of issues and possible generalizations to be resolved in 
concrete problems. 

One important generalization is a physically most interesting limit of a vanishing bulk mass 
M 2 , whose rigorous treatment should justify a qualitative comparison of the above type for the 
cutoff scales in our expansion (|72p and the weak field expansion of the DGP model. The local 
curvature expansion is perfect and nonsingular for nonvanishing M 2 and is applicable within 
its cutoff scale (fl"6|) . However, for M 2 — > it obviously breaks down, because the proper time 
integrals start diverging at the upper limit and all UV finite terms of (|15p blow up. These infrared 
divergences can be avoided by a nonlocal curvature expansion of the heat kernel of [23]. Up to 
the cubic order in curvatures this expansion explicitly exists for Tr e s ^ |24j , but for the structure 
involving a local differential operator TrVF(V)e sD it still has to be developed. 

Another important generalization is the extension of these calculations to the cases when 
already the lowest order approximation involves a curved spacetime background (i.e. dS or AdS 
bulk geometry, deSitter rather than flat brane, etc.). The success of the above technique is 
obviously based on the exact knowledge of the y-dependence in the lowest order Green's function 
in the bulk and the possibility to perform exactly (or asymptotically for large M 2 ) the integration 
over y. All these generalizations and open issues are currently under study. 

To summarize, we developed a new scheme of calculating quantum effective action for the 
braneworld DGP-type system in curved spacetime. This scheme gives a systematic curvature 
expansion by means of a manifestly covariant technique. Combined with the method of fixing 
the background covariant gauge for diffeomorphism invariance developed in \12\ [25] this gives the 
universal background field method of the Schwinger-DeWitt type in gravitational brane systems. 
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